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ABSTRACT 

We examine what types of galaxies the conversion formula from dust infrared (IR) luminosity into the 
star formation rate (SFR) derived by Kennicutt (1998) is applicable to. The ratio of the observed IR 
luminosity, £ir, to the intrinsic bolometric luminosity of the newly (< 10 Myr) formed stars, £sf, of 
a galaxy can be determined by a mean dust opacity in the interstellar medium and the activity of the 
current star formation. We find that these parameters area being 0.5 < Lir/Lsf < 2.0 is very large, and 
many nearby normal and active star-forming galaxies really fall in this area. It results from offsetting 
two effects of a small dust opacity and a large cirrus contribution of normal galaxies relative to starburst 
galaxies on the conversion of the stellar emission into the dust IR emission. In conclusion, the SFR 
determined from the IR luminosity under the assumption of Lir = L$f like Kennicutt (1998) is reliable 
within a factor of 2 for all galaxies except for dust rich but quiescent galaxies and extremely dust poor 
galaxies. 

Subject headings: dust, extinction — galaxies: general — infrared: galaxies — stars: formation 



1. INTRODUCTION 

Dust infrared (IR) luminosity 1 is a familiar indicator 
of the current star formation rate (SFR). Especially, the 
most active star- forming galaxies (e.g., Arp 220) emit al- 
most all their radiation energy in the IR range (Soifer, 
Houck, & Neugebauer 1987; Sanders & Mirabel 1996). For 
such galaxies, Kennicutt (1998a) derived a conversion for- 
mula from the IR luminosity to the SFR by assuming the 
IR luminosity equal to the bolometric luminosity and 10- 
100 Myr continuous star formation (hereafter K98 IR for- 
mula). More recently, Inoue, Hirashita, & Kamaya (2000; 
hereafter IHKOO) theoretically constructed a new conver- 
sion law from IR to SFR, which we can apply not only to 
dusty starbursts but also to normal (or quiescent) galaxies. 

Although the SFRs determined from various indicators 
should agree with each other, we often find discrepancies 
among them. Thus, many authors try to resolve such dis- 
crepancies (Sullivan et al. 2000; Bell & Kennicutt 2001; 
Hopkins et al. 2001; Chariot et al. 2002; Buat et al. 2002; 
Rosa-Gonzalez et al. 2002). In this context, the SFR from 
the IR luminosity via K98 IR formula is assumed to be the 
correct SFR in spite of the fact that K98 IR formula can be 
applied to only dusty starbursts in principle (Hopkins et 
al. 2001; Rosa-Gonzalez et al. 2002). Indeed, the SFR via 
K98 IR formula shows a good agreement with that from 
the Ha luminosity even for samples of normal galaxies if 
we correct the Ha luminosity for the dust extinction in 
the interstellar medium (ISM), that in H n regions, and 
the stellar Balmer absorption (Chariot et al. 2002; Rosa- 
Gonzalez et al. 2002). Also, an empirical conversion factor 
into the SFR from the IR luminosity for normal late-type 
spiral galaxies agrees with that of K98 IR formula within 
a factor of 2 (Buat & Xu 1996). Why can K98 IR formula 
present the correct SFR even for normal galaxies? 

In this Letter, we examine what types of galaxies K98 



IR formula is applicable to, and how much correction fac- 
tor against the formula is required to obtain a more precise 
SFR. In §2, we formulate the relation between the IR lu- 
minosity and the star formation (SF) luminosity. In §3, 
we describe a model of dust extinction of galaxies. In §4, 
we show that the ratio of the IR luminosity to the SF lu- 
minosity becomes of order unity within a factor of 2 for a 
very wide parameter area, and many nearby galaxies really 
fall in this area. 

2. INFRARED LUMINOSITY OF GALAXIES 

Following IHKOO, we formulate the relation between the 
observed IR luminosity of a galaxy, Lir, and the intrin- 
sic bolometric luminosity of the newly (< 10 Myr) formed 
stellar population (hereafter SF population), Lsf- 

The IR luminosity of a galaxy is expressed by the fol- 
lowing: 



Lm = L Lya + (1 - f)L LC + e^uv + f]L 



under ? 



(i) 



where LLya, Llc, £uv> and L U nder are luminosities of Ly- 
man a photons generated in H n regions, Lyman contin- 
uum (LC) photons and nonionizing (A > 912 A) photons 
emitted by the SF population, and photons from underly- 
ing stellar population, respectively. The parameters, /, e, 
and rj are a fraction of LC photons contributing to hydro- 
gen ionization, an efficiency of the dust extinction for non- 
ionizing photons from the SF population, and that for pho- 
tons from underlying population. Here, we assume that all 
Lyman a photons are eventually absorbed by dust during 
a number of resonance scatterings in the ISM. We neglect 
the escape of LC photons from a galaxy, the consumption 
of LC photons by Helium ionization, any extinctions for 
hydrogen's recombination lines within H n regions except 
for Lyman a, and the contribution to the IR luminosity 
by AGN activity. 

In this Letter, the term 'IR' means the whole wavelength range of dust emission. The 'IR luminosity' dose not include the component emitted 
by stars directly. 
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According to the Kurucz ATLAS 9 stellar spectra with 
Solar metallicity and a turbulence speed of 2 km s^ 1 , we 
obtain Llc = 0.37Lsf and L\jy = 0.63Lsf if the Salpeter 
IMF (0.1-100 M Q ) is assumed. Since about two-thirds 
of the ionization-recombination processes produce Lyman 
a photons under case B (Spitzcr 1978), we also obtain 
Lhya = 0.12/Lsf- Here, we note that only a fraction 
/ of LC photons ionize hydrogen. Moreover, we intro- 
duce a new parameter, 7, to represent the intrinsic stel- 
lar luminosity fraction of the SF population to all stellar 
population; 7 = L SF /L totah where L to tai = -^sf + Sunder- 
This parameter is an indicator of the star-forming activity. 
Therefore, equation (1) is reduced to 



— = 0.37 - 0.25/ + 0.63e + tj— 

Lsf 7 



7 



(2) 



For dusty starbursts (7 ~ 1), the right hand side of equa- 
tion (2) becomes unity since the dust opacity is enough 
high to absorb all photons emitted by stars (i.e. / = 
and e = 1). In principle, K98 IR formula can be applied 
to only such galaxies. On the other hand, if we deter- 
mine all four parameters (/, e, 77, and 7), we can convert 
the observed IR luminosity of all kinds of galaxies to the 
SF luminosity, and then, to the current SFR. This is the 
algorithm of IHK00. 

3. DUST OPACITY OF GALAXIES 

Here, we describe a model of dust extinction to deter- 
mine the values of /, e, and 77, in equation (2). We will 
discuss in §4.2 the effect of the choice of a specific extinc- 
tion law on our result. 

Following Chariot & Fall (2000) (especially their Fig. 
1), we consider, first, two types of dust extinction. One is 
the extinction by diffuse ISM dust (r^ SM ), and the other 
is by dust accompanying the birth clouds of young stars 



BG\ 



Then, we assume rf M = Ty SM (A/Ay)~" and 



if c = nTy M (\/\ v )~ n for A > 912 A, where fj, = 2, 
n = 0.7, and Ay = 5500 A are adopted, so as to reproduce 
the observed relation between UV flux slope and lumi- 
nosity ratio of IR to UV for nearby star-forming galaxies 
(Chariot & Fall 2000). We define the dust optical depth 



BC 



.ISM 



for the SF population as t^ f = t- 

For simplicity, we assume that the dust geometry is 
an uniform screen (or the optical depth defined above 
is an effective one). Then, we define e as an average 
of e(i) over the age range of - 10 Myr, where e(t) = 

1 - J L x (t)cxp(-Tf F )dX/ J L x (t)d\ and L x (t) is the lu- 
minosity density of a simple stellar population (SSP) with 
an age t. We use SSPs produced by PEGASE 2.0 (Fioc 
& Rocca-Volmerange 1997) throughout this Letter. In the 
same way, we define r\ as an average of rj(t) over the age 
of 10 Myr - 15 Gyr but rf F is replaced with t{ sm . The 
obtained e and r\ are a function of only Ty SM , the normal- 
ization of the adopted extinction law. 

Let us introduce one more type of dust extinction to 
determine the parameter /, Lyman continuum extinction 
(Inoue, Hirashita, & Kamaya 2001, Inoue 2001), which is 
the extinction for LC photons by dust in H 11 regions be- 
fore LC photons ionize neutral hydrogens. The parameter 

2 The mean value is different from that in Inoue et al. (2001) because 
density at 200 (im (-30%; Alton ct al. 1998). In addition, we should 
extinction properties for a number of line of sights, which is different 



/ is determined from the dust optical depth for LC pho- 
tons in H 11 regions, r™, via eq. (8) in Petrosian, Silk, & 
Field (1972). Here, we assume a scaling law between t^q 



and 



.ISM 



i.e. 



r HII 

r LC 



/r- 



ISM 

V 



Now, we attempt to estimate £ from real data. We have 
shown that r™ is proportional to the dust-to-gas mass ra- 
tio, X>, in Inoue et al. (2001). Here, we adopt the relation, 
r^g ~ (D/Dmw) (Inoue 2001), where V M w = 6 x 10~ 3 
is a typical Galactic value (Spitzer 1978). We estimated 
T> for seven nearby spiral galaxies observed by ISO (Al- 
ton et al. 1998) in Inoue et al. (2001). The mean value is 
1.8 x 10~ 3 . 2 On the other hand, we can estimate the dust 
optical depth from the IR surface brightness. By neglect- 
ing the self-absorption of the IR photons and assuming 
the obtained optical depth to correspond to that in the 
diffuse ISM, we obtain t^ sm = 0.36 as a mean dust optical 
depth for these seven galaxies, where we adopt the dust 
emissivity determined by Bianchi, Davies, & Alton (1999). 
Therefore, we obtain £ « 1. 

Although the real relation between r™ and Ty SM is 
quite uncertain, we find in equation (2) that the effect of 
/ on the ratio Ltr/Lsf is small. Indeed, the effect of the 
deviation from / = 0.5 on the Lir/Lsf is at most ±0.125. 

4. LUMINOSITY RATIO OF DUST INFRARED EMISSION TO 
STAR-FORMING POPULATION 

Now, we can determine LmJ L$f from only two pa- 
rameters; a mean dust opacity (ry' M ) and a star-forming 
activity (7). We survey these parameters space to clar- 
ify when we can determine the SFR from Lir within a 
factor of 2 uncertainty via K98 IR formula, i.e., when 
0.5 < Lir/Lsf < 2.0 is satisfied. 

In Figure 1, we show some constant Lir/Xsf lines in 
the Ty SM -7 plane. We define the 'applicable area' as the 
parameter area satisfying the condition 0.5 < Lir/Lsf < 
2.0. The 'applicable area' is between two thick lines in Fig- 
ure 1 since the upper thick line is the case of Lir/Lsf = 0.5 
and the lower one is the case of Lir/Lsf = 2.0. Also, the 
thin line in Figure 1 represents the case of Lir/Lsf = 1-0. 
We find that the 'applicable area' is very wide in the Ty SM - 
7 plane. This is a main conclusion in this Letter. 

For a galaxy with an exponentially declining star for- 
mation history (SFH) of > 10 Gyr time scale, age of > 1 
Gyr, and Salpeter IMF, we find that 7 ~ 0.5. Thus, many 
spiral and irregular galaxies are likely to fall in the 'appli- 
cable area'. In order to examine whether real galaxies fall 
in the 'applicable area' or not, we will estimate their 7 and 
t ism m j. ne following. Then, we will discuss implications 
from our results. 

4.1. Estimation 0/7 and Ty M 

First, we estimate 7 of galaxies from the observed equiv- 
alent width of H a emission line corrected for dust and 
stellar absorption effects. By the definition of 7, we can 
express 7 = EWn a {Cn a / C r) where EWna is the intrinsic 
equivalent width of H a emission line, Ch q is the intrinsic 
ratio of the bolometric luminosity to the H a emission line 
luminosity of the SF population, and Cr is the intrinsic ra- 
tio of the bolometric luminosity to the luminosity density 

we correct it for He contribution and for the calibration error of the flux 
note that X>mw by Spitzer (1978) is estimated from the analysis of the 
from our method. 



SFR from IR 



3 



at i?-band of all stellar population. 

If we specify the metallicity and IMF of galaxies, the 
value of CW is determined from only the electron tem- 
perature and density in H II regions under case B approx- 
imation. On the contrary, the value of Cr depends on 
the SFH and age of galaxies as well as the metallicity and 
IMF. We examine the evolution of Cr with the galactic 
age, so that we find that the factor Chq/Cr falls in the 
range of 0.002 - 0.004 unless we consider a galaxy younger 
than a few Gyr or a galaxy with a short (< 1 Gyr) ex- 
ponential time scale. Thus, we adopt a constant 0.003 as 
Chq/Ck, which corresponds to the case of Solar metallic- 
ity, exponentially declining SFH with 10 Gyr time scale, 
and about 13 Gyr age. Here, we do not take into account 
the chemical enrichment effect. 

Therefore, we estimate 7 from the following equation 
approximately; 

f 0.003£W Ha for EW Ha < 333 A ,^ 
7 \ 1 for EW na > 333 A ' [ ' 

where the second case is needed due to 7 < 1 by its def- 
inition. In the real case, the value of 7 asymptotically 
approaches unity by decreasing the factor Ch q /Cr, when 
EWn a increases. Thus, we overestimate 7 for a galaxy 
with a very large EWn a ■ However, the real value of 7 is 
about unity for such galaxies. On the other hand, we un- 
derestimate 7 for a galaxy with a very small EWn a . We 
confirm that this effect is also not so large (at most factor 
of 2). 

Although the effect of the change of the metallicity 
on Cua/Cn is small, that of the IMF slope is large. 
If we choose the IMF slope of 1.35 (Salpeter's slope is 
2.35), a typical Cn a /Cn decreases an order of magnitude, 
about 0.0003. In the case of the slope of 3.35, we find 
Cua/Cfs, — 0.007 typically. This is because a top-heavy 
IMF (1.35 slope) results in a smaller fraction of the R- 
band luminosity density in the bolometric luminosity since 
massive stars do not contribute to the i?-band light so ef- 
fectively. A steep IMF provides an opposite result. 

Now, we must estimate the intrinsic equivalent width, 
EWna, from the observed one, EW^. Assuming that the 
underlying stellar population (age > 10 Myr) dominates 
the i?-band continuum and the stellar Balmer absorption, 
we can express 

EW Ua = {EW$£ + EWZ hs )j exp(r A s H F a - r^ 1 ) , (4) 

where EW^ hs is the equivalent width of the stellar ab- 
sorption, and Ah q = 6563 A. The factor 1// means that 
only a fraction / of LC photons is used to ionize neutral 
hydrogen atoms. 

Next, we estimate t v f from the Balmer decrement cor- 
rected for the stellar absorption by assuming EW^ hs = 
EW£p hs = 5 A, which noted by Kcnnicutt (1992) as a 
typical value, and then, we convert the obtained t v f into 
Ty SM by using the parameter /x = 2 (See §3). Here, we 
also assume the electron temperature of 10 4 K and case 
B. Once the optical depth Ty M is obtained, we can esti- 
mate / parameter, EW^a from EWffa, and then, 7. The 
obtained values of Ty M are systematically smaller than 
but roughly consistent with those determined from the 
UV slope method (Meurer, Heckman, & Calzctti 1999) 
for sample galaxies described below. 



The uncertainty of Ty from the Balmer decrement is 
large, ~ ±0.2. As a result, the uncertainty of the obtained 
7 is <~ ±0.2, which is estimated from a simple Monte Carlo 
simulation if we do not take account of the uncertainties 
of the adopted parameters of the dust extinction law (i.e., 
77, £, and n in §3). Also, we find that a random error 
of Ty M causes a somewhat overestimation of 7. To dis- 
play this point, we show a typical error bar in Figure 1. 
In addition, we select only galaxies with enough strong H 
(3 'emission' line to obtain a reliable Balmer decrement. 
Thus, our sample galaxies are biased towards larger 7. 

One might think that the parameter 7 can be estimated 
from the ratio of the UV luminosity and the R or iVband 
luminosity. However, we can do it reasonably only if we 
specify the SFH and age of galaxies. That is, the coefficient 
between 7 and Ljjv /Lr or Ljjv/Lk varies together with 
the galactic age significantly. This is because a significant 
fraction of the UV light originates in the underlying stellar 
population with its age of > 10 Myr. On the contrary, the 
ionizing photons producing the H a photons originate in 
almost only the SF population (its age < 10 Myr). There- 
fore, we estimated 7 from the equivalent width of the H a 
line. 



4.2. Results & Discussions 

To plot real galaxies in Figure 1 , we compile all morpho- 
logical types of galaxies except for elliptical galaxies from 
the literature. In Figure 1, the filled circles are early type 
spiral (Sa-Sab) galaxies observed by Usui, Saito, & Tomita 
(2001), who selected active star-forming early type spi- 
ral galaxies with log(LFiR/is) > 0.5. The filled squares, 
open squares, and crosses are early type spiral (SO-Sbc), 
late type spiral (Sc-Sm), and irregular (including peculiar 
and merger) galaxies, respectively, observed by Kennicutt 
(1992), who selected galaxies covering all morphological 
types. The filled triangles, open triangles, and stars are 
also early spiral, late spiral, and irregular (also peculiar 
and merger) galaxies, respectively, observed by Storchi- 
Bergmann, Kinney, & Challis (1995), who selected ac- 
tive star-forming galaxies observed by IUE satellite. Any 
galaxies containing an AGN are excluded. 

We find that almost all (43/52) galaxies fall in the 'ap- 
plicable area'. For 7 out of 9 galaxies in the outside of the 
'applicable area', we fail to determine Ty M because they 
become negative values. Although the uncertainties of the 
determined Ty M and 7 are large and compiled galaxies 
do not make a complete sample, this result indicate that 
many nearby normal galaxies are really in the 'applicable 
area'. 

We emphasize the following two points. First, it is a just 
coincidence that K98 IR formula is applicable to normal 
galaxies. It results from offsetting two effects of a small 
dust opacity and a large cirrus contribution of normal 
galaxies relative to starburst galaxies on the conversion 
coefficient of K98 IR formula (Kennicutt 1998b, IHK00). 
For an example case of a galaxy with Ty M = 0.3 and 
EW^ = 50 A, while only about 70% luminosity of the 
SF population is converted into the IR luminosity, it is 
almost the same amount as the rest 30% SF luminosity 
that is supplied with the IR luminosity by the underlying 
population. 
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Second, the result above is robust for changing the 
adopted parameter sets of the dust extinction. One might 
think that our result depends on the dust extinction law 
adopted. However, we note that it dose not change when 
we change the parameters in the range of < /i < 3, 
0.4 < n < 1.0, and < £ < 2 as shown in Figure 2 (an 
extra figure). This point will be discussed in detail in the 
forthcoming paper (Inoue 2002, in preparation). 3 

We discuss the condition required to fall in the 'ap- 
plicable area'. From the lower bound of 7 in Figure 1 
(Lir/Lsf = 2.0), we find that a quiescent galaxy with a 
large dust amount (a small 7 and a large Ty M ) is in the 
outside of the 'applicable area'. Since galaxies showing a 
lower activity of the star formation tend to be more trans- 
parent (Wang & Heckman 1996; Hopkins et al. 2001), the 
number of galaxies falling in the area of large Ty M and 
small 7 may be small. Moreover, we find that the lower 
bound of 7 is ~ 0.1 for a galaxy with r^ SM - 0.3. This 
corresponds to EW^ > 10 A. This boundary can be ex- 
ceeded not only by active star-forming galaxies but also by 
many nearby normal spiral and irregular galaxies (Kenni- 
cutt 1983). 

The upper bound of 7 in Figure 1 (Lir/£sf = 0.5) 
reaches unity if Ty is larger than about 0.1. Thus, a star- 
burst galaxy having very small dust content of M < 0.1 
is in the outside of the 'applicable area'. For example, I Zw 
18, which is an active star-forming and the most metal de- 
ficient nearby dwarf galaxies, has typically Ay = 0.2 mag 



in the southeast H 11 region (Cannon et al. 2002). It cor- 
responds to Ty M ~ 0.06 if we adopt fi = 2. By taking into 
account 7 ~ 1 estimated from very large EW^ ~ 1000 
A (Cannon et al. 2002), I Zw 18 is in the outside of the 
'applicable area'. 

In conclusion, we find that the range of a dust opacity 
and a star- forming activity to be 0.5 < Ltr/£sf < 2.0 is 
very large. Furthermore, we confirm that these parame- 
ters of many nearby normal and active star-forming galax- 
ies really fall in this range. Therefore, we can apply K98 
IR formula assuming Lir = L$f to many nearby normal 
galaxies as well as dusty starbursts without any correc- 
tion factor within a factor of 2 uncertainty. However, we 
cannot apply it to dust rich but quiescent galaxies and ex- 
tremely dust poor galaxies. Therefore, we need to divide 
the coefficient of K98 IR formula by a correction factor, 
Lir/Lsf, when we apply it to such galaxies. Following the 
model in this Letter, we can estimate the correction factor 
for any galaxies if we have all data required, although we 
must more refine the model. 
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Fig. 1. — Iso-Lir/Lsf lines 



. the - 



-7 plane, where Ltr is the observed total IR luminosity, Lgp is the intrinsic bolomctric luminosity of 



newly (< 10 Myr) formed stellar population, tJP m is a mean optical depth by dust in the diffuse interstellar medium, and 7 is the luminosity 
fraction of newly formed stellar population. The upper thick line is Lj^/Lgp = 0.5, and the lower thick line is Lj^/Lgp = 2.0. The thin 
line is Ljj^/Lgp = 1.0. The filled circles are early type spiral (Sa-Sab) galaxies taken from Usui, SaitS, & Tomita (2001). The filled squares, 
open squares, and crosses are early type spiral (SO-Sbc), late type spiral (Sc-Sm), and irregular (including peculiar and merger) galaxies, 
respectively, taken from Kcnnicutt (1992). The filled triangles, open triangles, and stars are also early, late spirals, and irregular galaxies, 
respectively, taken from Storchi-Bcrgmann, Kinney, & Challis (1995). The classifications arc taken from NASA/IPAC Extragalactic Database. 
No galaxies with an AGN arc included. A typical error bar is shown in the bottom right of the panel. 
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